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1. INTRODUCTION 

The term "QCD vacuum" is frequently abused. 
Only in the case of the Hamiltonian formulation 
is it clear what we mean by the vacuum: it is the 
wave functional associated with the lowest energy 
state. Observables create excitations on top of 
this vacuum. Knowing the vacuum is knowing 
all: We should know better. 

Strictly speaking the vacuum is empty. Nev- 
ertheless its wave functional can be highly non- 
trivial, deviating considerably from that of a non- 
interacting Fock space, based on a quadratic the- 
ory. Even in the later case the result of probing 
the vacuum by boundaries is non-trivial as we 
know from Casimir. The probe is essential: one 
needs to disturb the vacuum to study its prop- 
erties. Somewhat perversely the vacuum may be 
seen as a relativistic aether. It promises to mag- 
ically resolve our problems, from confinement to 
the cosmological constant. For the latter super- 
symmetry is often called for to remove the other- 
wise required fine-tuning. It merely hides the rel- 
ativistic aether, even giving it further structure. 
Remarkably it seems to have enough structure to 
give a non-trivial example of the dual supercon- 
ductor at work [Q. 

Most will indeed put their bet on the dual su- 
perconductor picture for the QCD vacuum [Q, 
and this has motivated the hunt for magnetic 
monopoles using lattice techniques, long before 
supersymmetric duality stole the show [Q. The 
definitions rely on choosing an abelian projec- 
tion 1^ and the evidence is based on the no- 
tion of abelian dominance Q, establishing the 
dual Meissner effect H, or the construction of 



a magnetically charged order parameter, whose 
non-zero expectation value implies spontaneous 
breaking of the dual gauge symmetry , yielding 
electric confinement. But center vortices, probed 
by the newly defined center dominance, this year 
suddenly became center stage again and we 
will surely hear more next year. 

In the euclidean path integral only the vacuum 
state will contribute when we let the (imaginary) 
time go to infinity, underlying the essence of the 
transfer matrix approach to extract observables 
from euclidean path integrals as used in lattice 
gauge theories. However, this is not what peo- 
ple have in mind when talking about the vac- 
uum structure of gauge theories. More appro- 
priate for most studies would be to talk about 
the low-energy and long-distance behaviour of the 
theory. One way to address this is to attempt to 
isolate the relevant degrees of freedom for which 
one can derive an effective low-energy theory. 
Monopole actions derived from block-spin trans- 
formations 1^ and the instanton liquid model ||] 
are examples. It is not required that the relevant 
degrees of freedom are associated to semi-classical 
objects, the main reason being that in a strongly 
interacting theory the quantum fluctuations can 
be (much) bigger than the classical effects. Some- 
times this argument is used against the relevance 
of instantons, which by their very definition as 
localised objects in time might lose significance 
when quantum fluctuations are too strong. 

Due to the limited space available this re- 
view concentrates on the instanton content of the 
theory, where recently considerable progress has 
been achieved. But first I will use this opportu- 
nity to explain my own thoughts on the matter. 
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2. VACUUM DEMOCRACY 

The model we wish to describe here starts from 
the physics in a small volume, where asymptotic 
freedom guarantees that perturbative results are 
valid. The assumption is made, that at least for 
low-energy observables, integrating out the high- 
energy degrees of freedom is well-defined pertur- 
batively and all the non-perturbative dynamics is 
due to a few low- lying modes. This is most easily 
defined in a Hamiltonian setting, since we are in- 
terested in situations where the non-perturbative 
effects are no longer described by semiclassical 
methods. 

2.1. Complete gauge fixing 

Due to the action of the gauge group on the 
vector fields, a finite dimensional slice through the 
physical configuration space (gauge inequivalent 
fields) is bounded. One way to demonstrate this 
is by using the complete Coulomb gauge fixing, 
achieved by minimising the norm of the gauge 
field along the gauge orbit. At small energies, 
fields are sufficiently smooth for this to be well 
defined and it can be shown that the space under 
consideration has a boundary, defined by points 
where the norm is degenerate. These are by defi- 
nition gauge equivalent such that the wave func- 
tionals are equal, possibly up to a phase factor in 
case the gauge transformation is homotopically 
non-trivial. The space thus obtained is called a 
fundamental domain. For a review see ref. |p^ . 

2.2. Non-perturbative dynamics 

Given a particular compact three dimensional 
manifold M on which the gauge theory is defined, 
scaling with a factor L allows one to go to larger 
volumes. It is most convenient to formulate the 
Hamiltonian in scale invariant fields A — LA. Di- 
viding energies by L recovers the L dependence in 
the classical case, but the need of an ultraviolet 
cutoff and the resulting scale anomaly introduces 
a running coupling constant g{L), which in the 
low-energy effective theory is the only remnant of 
the breaking of scale invariance. 

When the volume is very small, the effective 
coupling is very small and the wave functional is 
highly localised, staying away from the bound- 
aries of the fundamental domain. We may com- 



pare with quantum mechanics on the circle, seen 
as an interval with identifications at its bound- 
ary. At which points we choose these boundaries 
is just a matter of (technical) convenience. The 
fact that the circle has non-trivial homotopy, al- 
lows one to introduce a 9 parameter (playing the 
role of a Bloch momentum) . 

Expressed in A, the shape of the fundamental 
domain and the nature of the boundary condi- 
tions, is independent of L. Due to the rise of the 
running coupling constant with increasing L the 
wave functional spreads out over the fundamental 
domain and will start to feel the boundary iden- 
tifications. This is the origin of non-perturbative 
dynamics in the low-energy sector of the theory. 

Quite remarkably, in all known examples (for 
the torus and sphere geometries), the sphalerons 
lie exactly at the boundary of the fundamental 
domain, with the sphalcron mapped into the anti- 
sphaleron by a homotopically non-trivial gauge 
transformation. The sphaleron is the saddle point 
at the top of the barrier reached along the tun- 
nelling path associated with the largest instanton, 
its size limited by the finite volume. 

For increasing volumes the wave functional first 
starts to feel the boundary identifications at these 
sphalerons, "biting its own tail". When the en- 
ergy of the state under consideration becomes of 
the order of the energy of this sphaleron, one can 
no longer use the semiclassical approximation to 
describe the transition over the barrier and it is 
only at this moment that the shift in energy be- 
comes appreciable and causes sizeable deviations 
from the perturbative result. This is in particular 
true for the groundstate energy. Excited states 
feel these boundary identifications at somewhat 
smaller volumes, but nodes in their wave func- 
tional near the sphalcron can reduce or postpone 
the influence of boundary identifications. 

This has been observed clearly for SU{2) on a 
sphere . The scalar and tensor glueball mass is 
reduced considerably due to the boundary identi- 
fications, whereas the oddball remains unaffected 
(see fig. 1). These non-perturbative effects re- 
move an unphysical near-degeneracy in perturba- 
tion theory (with the pseudoscalar even slightly 
lower than the scalar glueball mass). The dom- 
inating configurations involved are associated to 
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instanton fields, in a situation where semiclassical 
techniques are inappropriate for computing the 
magnitude of the effect. When boundary identi- 
fications matter, the path integral receives large 
contributions from configurations that have non- 
zero topological charge, and in whose background 
the fermions have a chiral zero mode, its conse- 
quences to be discussed later. 





mR 

:t 














11 


> f 

1 1 




0.1 


0.25 



Figure 1. The low- lying glueball spectrum on a 
sphere of radius i? as a function of / = g^(i?)/27r^ 
at 9 = 0. Approximately, / = 0.28 corresponds to 
a circumference of 1.3 fni. From ref. [11]. 

At some point technical control is lost, since 
so far only the appropriate boundary conditions 
near the sphalerons can be implemented. As soon 
as the wave functional starts to become apprecia- 
ble near the rest of the boundary too, this is no 
longer sufficient. 

This method has in particular been very suc- 
cessful to determine the low-lying spectrum on 
the torus in intermediate volumes, where for 
SU{2) agreement with the lattice Monte Carlo re- 
sults has been achieved within the 2% statistical 
errors jl^Jl^. In this case the non-perturbative 
sector of the theory was dominated by the en- 
ergy of electric fiux (torelon mass), which van- 
ishes to all orders in perturbation theory. The 
leading semiclassical result is exp(— S'o/.g(i)), due 
to tunnelling through a quantum induced barrier 
of height £'s = 3.21 /i and action S'o = 12.5. Al- 
ready beyond 0.1 fm this approximation breaks 
down. One finds, accidentally in these small vol- 
umes, the energy to be nearly linear in L. 

The effective Hamiltonian in the zero-momen- 
tum gauge fields, derived by Liischer |l^, and 
later augmented by boundary identifications to 
include the non-perturbative effects llj], breaks 



down at the point where boundary identifica- 
tions in the non-zero momentum directions as- 
sociated with instantons become relevant. The 
sphaleron has an energy 72.605/(Lg^(L)) and 
was constructed numerically |lj]. Its effect be- 
comes noticeable beyond volumes of approxi- 
mately (0.75 fm)3. For SU{3) this was verified 
directly in a lattice Monte Carlo calculation of 
the finite volume topological susceptibility [ |T5| . 
The results for the sphere have shown that also 
these effects can in principle be included reliably, 
but the lack of an analytic instanton solution on 
T'^ X IR has prevented us from doing so in practise. 

2.3. Domain formation 

The shape of the fundamental domain depends 
on the geometry. Assuming that g{L) keeps on 
growing with increasing L, causing the wave func- 
tional to feel more and more of the boundary, one 
would naturally predict that the infinite volume 
limit depends on the geometry. This is clearly 
unacceptable, but can be avoided if the ground 
state obtained by adiabatically increasing L is not 
stable. Thus we conjecture that the vacuum is 
unstable against domain formation. This is the 
minimal scenario to make sure that at large vol- 
umes, the spectrum is independent of its geome- 
try. Domains would naturally explain why a non- 
perturbative physical length scale is generated in 
QCD, beyond which the coupling constant will 
stop running. However, we have no guess for the 
order parameter, let alone an effective theory de- 
scribing excitations at distances beyond these do- 
mains. Postulating their existence, nevertheless a 
number of interesting conclusions can be drawn. 

The best geometry to study domain formation 
is that of a box since it is space-filling. We can 
exactly fill a larger box by smaller ones. This is 
not true for most other geometries. In small to 
intermediate volumes the vacuum energy density 
is a decreasing function |l^ of L, but in analogy 
to the double well problem one may expect that 
at stronger coupling the vacuum energy density 
rises again with a minimum at some value Lq, 
assumed to be 0.75 fm. For L sufficiently larger 
than Lq it thus becomes energetically favourable 
to split the volume in domains of size Lq . 

Since the ratio of the string tension to the 
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scalar glueball mass squared shows no structure 
around (0.75 fm)^, we may assume that both have 
reached their large volume value within a domain. 
The nature of their finite size corrections is suffi- 
ciently different to expect these not to cancel ac- 
cidentally. The colour electric string arises from 
the fact that flux that enters the box has to leave 
it in the opposite direction. Flux conservation 
with these building blocks automatically leads to 
a string picture, with a string tension as com- 
puted within a single domain and a transverse 
size of the string equal to the average size of a 
domain, 0.75 fm. The tensor glueball in an in- 
termediate volume is heavily split between the 
doublet {E~^) and triplet {T^) representations of 
the cubic group, with resp. 0.9 and 1.7 times the 
scalar glueball mass. This implies that the tensor 
glueball is at least as large as the average size of 
a domain. Rotational invariance in a domain-like 
vacuum comes about by averaging over all orien- 
tations of the domains. This is expected to lead to 
a mass which is the multiplicity weighted average 
of the doublet and triplet, yielding a mass of 1.4 
times the scalar glueball mass. Domain formation 
in this picture is driven by the large field dynam- 
ics associated with sphalerons. Which state gets 
affected most depends in an intricate way on the 
behaviour of the wave functionals (cmp. fig. 1). 

In the four dimensional euclidcan context, 0(4) 
invariance makes us assume that domain forma- 
tion extends in all four directions. As is implied 
by averaging over orientations, domains will not 
neatly stack. There will be dislocations which 
most naturally are gauge dislocations. A point- 
like gauge dislocation in four dimensions is an in- 
stanton, lines give rise to monopoles and surfaces 
to vortices. In the latter two cases most natu- 
rally of the Zn type. We estimate the density of 
these objects to be one per average domain size. 
We thus predict an instanton density of 3.2 fm^^, 
with an average size of 1/3 fm. For monopoles we 
predict a density of 2.4 fm^"^. 

If an effective colour scalar field will play the 
role of a Higgs field, abelian projected monopoles 
will appear. It can be shown |l6| that a monopole 
(or rather dyon) loop, with its U{1) phase rotat- 
ing Q times along the loop (generating an elec- 
tric field), gives rise to a topological charge Q. 



In abelian projection it has been found that an 
instanton always contains a dyon loop ||l^. We 
thus argue this result to be more general, leading 
to further ties between monopoles and instantons. 

2.4. Regularisation and 6 

It is useful to point out that the non-trivial ho- 
motopy of the physical configuration space, like 
non-contractable loops associated to the instan- 
tons {Tri{A/G) = T^siG) = 7L), is typically de- 
stroyed by the regularisation of the theory. This is 
best illustrated by the example of quantum me- 
chanics on the circle. Suppose we replace it by 
an annulus. As long as the annulus does not fill 
the hole, or we force the wave function to vanish 
in the middle, theta is a well-defined parameter 
associated to a multivalued wave function. We 
may imagine the behaviour for small instantons 
in gauge theories to be similar to that at the cen- 
ter in the above model. Indeed, the gauge in- 
variant geodesic length of the tunnelling path for 
instantons on M x IR, given hy £ = dty^V{t), 
where V{t) is the classical potential along the tun- 
nelling path, is expected to vanish for instantons 
that shrink to zero size. This is confirmed for 
M = S^, using results of ref. ||l^ to show that 
(for unit radius) £^ 2r2(f)^67r/fe(l + 0(fo-i)), 
with the instanton size defined by p=(l + b^)~2 . 

Due to the need of regularising the ultraviolet 
behaviour, the small instantons are cut out of the 
theory. Using the lattice regularisation this does 
not leave a hole, but rather removes the "sin- 
gularity at the origin" , as the lattice configura- 
tion space has no non-contractable loops. Strictly 
speaking this means one can not have a theta pa- 
rameter at any finite lattice spacing. Furthermore 
the regularisation forces one to divide out all the 
gauge transformations as there are no homotopi- 
cally non-trivial ones. It is advisable to divide 
out all gauge transformations, even if some of the 
homotopy is preserved by some regularisation! 

One can, however, still introduce the theta pa- 
rameter by adding iOQ (with Q = J d^xq^x), 
q{x) = Tr{F^"'{x)Ff,^{x))/16TT'^ the topological 
charge operator) to the action. Of course only 
for smooth fields the charge Q will be (approx- 
imately) integer. Also within the Hamiltonian 
formulation one may introduce a theta parame- 
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ter by mixing the electric field with 6 times the 
magnetic field, E ^ E — 9B/{2ny. In these ap- 
proaches theta is simply a parameter added to the 
theory. Whether or not one will retrieve the ex- 
pected periodic behaviour in the continuum limit 
becomes a dynamical question. 

It should be pointed out that in particular for 
SU{N) gauge theories in a box (in sectors with 
non-trivial magnetic fiux) there is room to argue 
for a 2TrN, as opposed to a 27r, periodicity for 
the 9 dependence. However, the spectrum is pe- 
riodic with a period 2n, and the apparent dis- 
crepancy is resolved by observing that there is a 
non-trivial spectral flow This may lead to 

phase transitions at some value(s) of 9, related 
to the oblique confinement mechanism Sim- 
ilarly for supersymmetric gauge theories this in- 
terpretation, supported by the recent discovery 
of domain walls between different vacua |2^, re- 
moves the need for semiclassical objects with a 
charge 1/N. Such solutions do exist for the torus, 
but the fractional charge is related to magnetic 
flux and the interpretation is necessarily as stated 
above! The "wrong" periodicity in theta has long 
been used to argue against the relevance of in- 
stantons, but in the more recent literature this is 
now phrased more cautiously [^l[|22|] . 

3. INSTANTONS 

Instantons are euclidean solutions responsible 
for the axial anomaly, breaking the C/a(1) sub- 
group of the U{Nf) X U{Nf) chiral symmetry for 
Nf flavours of massless fermions as dictated 
by the anomaly, 9^^'/7^75^'/(a;) = 2A^/(7(x). 
The breaking of Ua{^) manifests itself in the 
semiclassical computations through the presence 
of fermion zero modes, with their number and 
chirality fixed by the topological charge, through 
the Atiyah-Singer index theorem [Q. Integra- 
tion over the fermion zero modes leads to the so- 
called 't Hooft vertex or effective interaction . 
The integration over the scale parameter of the 
instanton ensemble is infrared dominated and a 
non-perturbative computation is desirable. 

In addition it is believed that the instantons are 
responsible for chiral symmetry breaking, where 
a chiral condensate is formed, which breaks the 



axial gauge group UA{Nf) completely. This spon- 
taneous breaking is dynamical and it is less well 
established that instantons are fully responsible. 
It is the basis of the instanton liquid model as 
developed by Shuryak over the years. For a com- 
prehensive recent review see ref. The details 
of the instanton ensemble play an important role. 
Only a hquid-like phase, as opposed to the di- 
lute or crystalline phases, will give rise to a chiral 
condensate. The model also makes a prediction 
for the average size and the topological suscep- 
tibility. In particular the latter quantity should 
be well-defined beyond a semiclassical approxima- 
tion. For large sizes the instanton distribution is 
exponentially cut-off and instantons do not give 
rise to an area law for the Wilson loop. When 
large instantons are more weakly suppressed the 
situation may differ , but a semiclassical anal- 
ysis in this case should not be trusted. 

Remarkably the topological susceptibility in 
pure gauge theories can be related to the rj' mass 
through the so-called Witten-Veneziano relation, 
fUml+my2ml)/2Nf = JdiX<T{q{x)qiO))>R 
= Xt, leading to the prediction xt ^ (180 MeV)''. 
This is based on the fact that the C/4(l) symme- 
try is restored in the planar limit [ p6|j27[ |, with 
Xt of order 1/N'^. From the requirement that in 
the presence of massless quarks Xt (and the theta 
dependence) disappears, the pure gauge suscep- 
tibility can be related to the quark-loop contri- 
butions in the pseudoscalar channel. Pole dom- 
inance requires the lightest pseudoscalar meson 
to have a mass squared of order 1/N. Relating 
the residue to the pion decay constant gives the 
desired result . The index R indicates the ne- 
cessity of equal-time regularisation ||2^ . A deriva- 
tion on the lattice using Wilson and staggered 
fermions was obtained in ref. |28| , making use 
of Ward-Takahashi identities. Finally, also the 
coarse grained partition function of the instan- 
ton liquid model allows one to directly deter- 
mine the 't Hooft effective Lagrangian ||2^, from 
which the Witten-Veneziano formula can be read 
off This formula is almost treated as the 

holy grail of instanton physics. It is important 
to realise that some approximations are involved, 
although it is gratifying there are three indepen- 
dent ways to obtain it ||2|,pi28l|2i . 
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3.1. Field theoretic method 

A direct computation of xt — Jd^x < q{x)q{0) > 
on the lattice requires a choice of discretisation 
for the charge density. A particularly simple 
one is |0| qUx) = -J2Tr{U^,{x)Uf,,ix))/16Tr\ 
where U^u{x) — ^e^ivpaUpaix) is the clover aver- 
aged plaquette P^v{x), formed by the four plaque- 
ttes that meet at x. Due to the short distance sin- 
gularities the operators require renormalisations. 
The lattice charge operator is conjectured to sat- 
isfy Ql = idx) = a4QZ(/3) + 0(a6), correct- 
ing for the fact that Ql need not be integer. In 
addition an additive renormalisation, associated 
to the contact term discussed before |2^, arises 
due to mixing of xl with operators with the same 
quantum numbers, xl = a^Xt2■^(/3)^-M(/?}K?(a'')■ 
To determine Z{0) and M{0) one takes a clas- 
sical configuration with a fixed topological charge 
Q |31 . Monte Carlo updates are rejected if they 



change the charge as determined by cooling and 
subsequent measurements oi Ql and xl allows 
one to extract Z and M . (In a sector with fixed 
charge Q we note that Xt = Q^/^j with V the vol- 
ume.) Over a certain range, independence of the 
starting configuration and the volume has been 
observed and was henceforth assumed Due 
to the need to fix Q to determine the renormali- 
sation factors, this is in a sense a hybrid method, 
and might also be sensitive to some of the prob- 
lems of cooling - to be discussed shortly. 

Considerable progress has been achieved, how- 
ever, by repeating the study with other choices 
for qL, based on smearing the links (iteratively 
replacing them by staples) . This considerably re- 
duces the value of M, greatly facilitating the ex- 
traction of a signal. On a 16^ lattice one finds 
Xt = (175(5) MeV)4 (for SU{i) aX (3 = 5.9, 6.0 
and 6.1) and Xt = (198(8) MeV)" (for SU{2) at 
/3 = 2.44, 2.5115 and 2.57). For the discussion on 
finite temperature, how the scale was set and for 
further details and references see ref. pl[ |. 

3.2. Cooling 

In the continuum, the Schwartz inequality im- 
plies that the action of any field configuration 
with charge Q is bounded by 87r^|(5|, and is 
reached by (anti-)selfdual solutions. In its sim- 
plest form, cooling aims at finding this min- 



imal action, using it to identify the topologi- 
cal charge. As all quantum fluctuations are re- 
moved, no renormalisations are required, such 
that Xt=<Q^> /V , where the average is over the 
Monte Carlo ensemble. Cooling can be achieved 
by putting /3 = oo in the standard Monte Carlo 
update (accepting updates only if they lower the 
action). The same result is obtained by a sort of 
congruent gradient method, which uses the lattice 
equations of motion to lower the action ||3^ . For 
SU{2) this method is deterministic and allows for 
estimating its rate of convergence . 

It is a remarkable and deep mathematical prop- 
erty of non-Abelian gauge theories on a compact 
four dimensional manifold that exact SU (N) so- 
lutions exist for any charge Q with 4A^|(5| pa- 
rameters (moduli) , in general describing position, 
size and colour orientation of |(3| pseudoparticles. 
Although a compact manifold breaks the scale 
invariance, generically instantons with arbitrary 
size exist, only limited by the finite volume. A 
notable exception is charge one instantons on T**. 
There is no problem in having smooth configura- 
tions of unit charge, but in an attempt to make 
them self-dual they shrink to a point . 

To understand its implications we take the time 
direction to infinity, in which case finite action 
forces tunnelling from vacuum to vacuum config- 
uration. The vacuum on a torus is, however, not 
unique. Periodic boundary conditions in the time 
direction force the vacua at both ends to be the 
same and as soon as one releases this constraint, 
exact vacuum to vacuum tunnelling solutions ex- 
ist. This was studied for 0(3) through the exact 
solutions on a cylinder |3^. For gauge theories 
it can be proven that twisted (for SU{2) "anti- 
periodic") boundary conditions [Tg| ] remove the 
obstruction, even at finite T. Large instantons 
have finite size effects H of 0{l/T). For suffi- 
ciently large volumes, assuming the instanton size 
will be cutoff dynamically, the effect is irrelevant, 
further helped by the fact that in a large volume 
almost all configurations have higher charges. 

As always, the continuum limit competes with 
the infinite volume limit. Even with presentday 
computer power, the remaining window is uncom- 
fortably small. The above finite size effect is usu- 
ally swamped by the cutoff effects. As the lattice 
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cutoff breaks the scale and rotational invariance, 
we would expect that the action is no longer con- 
stant on the continuum moduli space. Indeed, for 
a smooth instanton the Wilson action behaves as 
Sw{p = pI a)^'&TT'^{l-p^'^ /b+0{p^^)) and causes 
the instanton to shrink, until it becomes of the 
size of the cutoff and falls through the lattice. 
Cooling will first remove high-frequency modes 
and one is left with a slow motion along the mod- 
uli space, giving rise to a plateau in the cooling 
history, used to identify the topological charge. 
One will miss instantons smaller than some fixed 
value Pc- Assuming asymptotic freedom, one eas- 
ily shows that the error vanishes in the contin- 
uum limit. Note that by construction, the cooling 
method never will associate charge to a disloca- 
tion with an action smaller than 9Qtt^ jWN , the 
entropic bound, which would spoil scaling p^ . 

For extracting the size distribution, cooling and 
under-relaxed (or slow) cooling |^ is problematic 
as the size clearly will depend on where along 
the plateau one analyses the data The size 
distribution can be made to scale properly only 
at the expense of carefully adjusts the number of 
cooling steps when going to different /3. 

One can avoid loosing instantons under cooling 
by modifying the action such that the scaling vio- 
lations change sign for example by adding a 
2x2 plaquette to the Wilson action. This so-called 
over-improved action has the property that in- 
stantons grow under cooling, until stopped by the 
finite volume. Consequently it would still muti- 
late the size distribution. This can be avoided by 
improving the action so as to minimise the scaling 
violations |40|. A particularly efficient choice is 
the so-called five-loop improved [bLi) lattice ac- 
tion: S^Li = Em,n Cm,n Y.x,tLU - Pm.ii,ny{x)) , 

where Pmfi,nv{x) is the m x n plaquette and 
ci,i = 1, C2,2 = -TK, ci^2 ci,3 = ^ and 

^3,3 = liw- Overimproved cooling |^ is defined 
by ci,i = (4 — e)/3 and C2,2 = (e — l)/48, which 
is 0{a^) improved for e = 0. The 5Li is im- 
proved to ©(a**) and allows for a tiny action bar- 
rier of less than one permill of the instanton ac- 
tion at Pc = 2.3. Instantons of smaller size will 
be rapidly lost under cooling, but larger ones will 
stay forever and practically not change their size, 
of course using twisted boundary conditions. For 



p > 4 the action is at most a factor 10~^ larger 
than its continuum value. After thus eliminat- 
ing the cutoff effects, the finite size effects in the 
charge one sector are clearly observed, see fig. 2. 
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Figure 2. Cooling history [40] for SU{2) charge 
one (and two) instanton with twisted and periodic 
boundary conditions. Squares represent 5Li and 
circles over- improved cooling (£ = —1). 

It is important to note that improved cool- 
ing Q is used as a diagnostic tool; the config- 
urations are still generated by the standard Wil- 
son action. One may of course also use improved 
actions for this purpose, but the 5Li action was 
simply not tuned dynamically. One difficulty in 
extracting the instanton distribution is that a 
typical ensemble will have both instantons and 
anti-instantons. These interact, although the ac- 
tion only slightly differs from Stt^ times the num- 
ber of pseudoparticles, provided they are not too 
close. Close instanton anti-instanton (I-A) pairs 
will annihilate. Cooling sufficiently long also re- 
moves more widely separated pairs. As the cool- 
ing makes the configuration smooth, the total 
charge can be measured reliably even in the pres- 
ence of these pairs. For SU(2), lattices of sizes 
12"* and 123x36 at /3 = 2.4 and 2.5 with periodic 
and twisted boundary conditions, as well as 24'* 
at P — 2.6 with twisted boundary conditions were 
used. At (3 = 2.5 a physical volume of 1.02 fm 
across gives sizable finite size corrections. One 
obtains Xt = (200(15) MeV)"* with good scaling 
properties which agrees with ref. |3l|. 

Extracting the size of an instanton is based 
on identifying the pseudoparticles. Using among 
other things the known profile of an isolated in- 
stanton, five different definitions for the size p 
were used, which are all well correlated The 
resulting size distribution is given in fig. 3, com- 
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billing results for /3 = 2.4 (averaged over the differ- 
ent lattice types and boundary conditions after 20 
cooling sweeps) and for /3 = 2.6 (after 50 cooling 
sweeps). The solid curve is a fit to the formula 
P{p) oc p^/3exp(-(p/w)P), with w = 0.47(9) fm 
and p = 3(l), which at small sizes coincides with 
the semiclassical result |2^ . The peak of this dis- 
tribution occurs at p — 0.43(5) fm. Under pro- 
longed cooling, up to 300 sweeps, I-A instanton 
annihilations and in particular finite size effects 
in the charge one sector do affect the distribu- 
tion somewhat, but not the average size, which 
therefore seems to be quite a robust result. 




0.2 0.4 0.6 0.8 1 

P(fm) 

Figure 3. SU{2) instanton size distribution for 
(3 = 2.4 (squares) and 2.6 (crosses) in a volume 
1.44 fm across at lattice spacings a = 0.12 and 
0.06 fm. The dotted and dashed lines represent 
the cutoff at apc for both lattices. From ref. [40]. 

It would be advantageous if one could come up 
with a definition for the size that is related to a 
physical quantity, since now the notion is based 
on the semiclassical picture. This is neverthe- 
less appropriate for the comparison with the in- 
stanton liquid. The relatively large value of the 
average size as compared to that of 1/3 fm pre- 
dicted by the instanton liquid is a point of 
worry, typically leading to stronger interactions 
that may lead to a crystal (without chiral sym- 
metry breaking), rather than a liquid. Neverthe- 
less, in ref. it has been tested that the pseu- 
doparticles are homogeneously distributed with a 
density of 2 — 3 fm""* and occupying nearly half 
the volume. This is the case when only close I-A 
pairs have annihilated and therefore depends on 



the amount of improved cooling. It does, how- 
ever, show that the pseudoparticles are relatively 
dense (more so than assumed in the instanton liq- 
uid j|]). The value of 3.2 fm"'' for the density, de- 
rived earlier in the context of the domain picture 
is quite realistic in the light of these results. 

3.3. Smoothing 

Another method to study instantons on the 
lattice is based on the classical fixed point ac- 
tions El| , defined through the saddle point equa- 
tion S^(y)=miinuy{S^P{U)+KT{U,V)). It is 
obtained as the weak coupling limit of the block- 
ing transformation with a positive definite ker- 
nel kT{U,V), which maps a lattice {U} to {V}, 
coarser by a factor two. Reconstructing the fine 
from the coarse lattice is called inverse blocking. 

It can be shown to map a solution of the lattice 
equations of motion to a solution on the finer lat- 
tice with the same action. Iterating the inverse 
blocking, the lattice can be made arbitrarily fine, 
thereby proving the absence of scaling violations 
to any power in the lattice spacing |^. This 
classically perfect action still looses instantons be- 
low a critical size, which is typically smaller than 
a lattice spacing. For solutions this most likely 
happens at the point where the continuum inter- 
polation of the lattice field is ambiguous, causing 
the integer geometric charge to jump. For 
rough configurations that are not solutions, in- 
verse blocking typically reduces the action by a 
factor 32 and makes it more smooth. The fixed 
point topological charge is defined as the limit- 
ing charge after repeated inverse blockings. This 
guarantees no charge will be associated to disloca- 
tions (of any action below the instanton action) . 

The classical fixed point action, although op- 
timised to be short range, still has an infinite 
number of terms and finding a suitable trunca- 
tion is a practical problem. From examples of 
parametrisations, the success of reducing scaling 
violations in quantities like the heavy quark po- 
tential (tested by restoring rotational invariance) 
is evident, for recent reviews see ref. ]45|] . In prac- 
tise only a limited number of inverse blockings is 
feasible and the fixed point topological charge has 
to rely on a rapid convergence. The closer one is 
able to construct the fixed point action the bet- 
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ter this convergence is expected to be. For two di- 
mensional non-linear sigma models sufficient con- 
trol was achieved to demonstrate that more than 
one inverse blocking did not appreciably change 



the topological charge |43 



In four dimensional gauge theories, both find- 
ing a manageable parametrisation and doing re- 
peated inverse blockings is a major effort. It goes 
without saying that if no good approximation to 
the fixed point action is used, one cannot rely on 
its powerful theoretical properties. Studies of in- 
stantons for SU{2) gauge theories were performed 
in ref . [Q . A 48 term approximation to the fixed 
point action was used to verify the theoretical 
properties. The geometric charge was measured 
after one inverse blocking and it was shown that 
for Q = 1 the instanton action was to within a 
few percent from the continuum action (slightly 
above it due to finite size effects), whereas for 
Q = the action was always lower. Subsequently 
a simplified eight parameter form was used on 
which the instanton action was somewhat poorly 
reproduced, but such that the Q = 1 boundary 
stayed above the entropic bound for the action. 

A value of Xt = (235(10) MeV)-* was quoted 
on a 8^ lattice with physical volumes of up to 
1.6 fm, taking full advantage of the fact that fixed 
point actions can be simulated at rather large lat- 
tice spacings. However, < > measured on the 
coarse lattice was up to a factor 4 larger than 
on the fine lattice (for the two dimensional study 
much closer agreement was seen [^|). Further 
inverse blocking to check stability of the charge 
measurement was not performed. 

The same eight parameter action was used in 
ref. but they did their simulations on the 

fine lattice and performed an operation called 
smoothing: first blocking and then inverse block- 
ing. They changed the proportionality factor k 
in the blocking kernel, requiring that the saddle 
point condition is satisfied for the blocked lat- 
tice. Due to the change of k the properties of 
the fixed point action that inspired these authors 
can unfortunately no longer be called upon as a 
justification. This smoothing satisfies the prop- 
erties of cooling (the action always decreases and 
stays fixed for a solution) and should probably 
by judged as such. (See for further comments be- 



low.) They consider their study exploratory and 
concentrate on finite temperature near and be- 
yond the deconfinement transition. 

In ref. |^7| the number of terms to parametrise 
the fixed point action was extended to four powers 
of resp. the plaquette, a six-link and an eight-link 
Wilson loop. The latter was required to improve 
on the properties for the classical solutions. They 
achieved pc = 0.94, still considerably smaller than 
for S'sLi, and reproduced the continuum instan- 
ton action to a few percent for p> pc- To increase 
the quality of the fit a constant was added, which 
should vanish in the continuum limit (as it drops 
out of the saddle point equation). Possible rami- 
fications of this at finite coupling are not yet suf- 
ficiently understood. 

After one inverse blocking insufficient smooth- 
ing is achieved to extract the pseudoparticle po- 
sitions and sizes and further inverse blocking was 
considered computationally too expensive. Like 
in ref. jio), they also introduced a smoothing cy- 
cle, but now by blocking the fine lattice back to 
the coarse one. Such a cycle would not change the 
action when the blocking is indeed to the same 
coarse lattice. However, there are 2'^ different 
coarse sublattices associated to a fine one and in 
ref. the smoothing cycle involved blocking to 
the coarse lattice shifted along the diagonal over 
one lattice spacing on the fine lattice. Unlike in 
ref. Gfl, the smoothing cycle will be repeated. 
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Figure 4. Example of the smoothing after 1 and 
9 cycles, shown on the fine lattice. From ref. [47]. 

Although the fixed point nature of the action 
guarantees it is close to a perfect classical action 
it needs to be demonstrated that it preserves the 
topological charge at sufficiently large scales. For 
cooling this is argued from the local nature of 
the updates, not affecting the long distance be- 
haviour. Improved cooling is in this sense less 
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local than ordinary cooling, and it might seem 
that due to the rather compact Wilson loops in- 
volved in the parametrisation of the fixed point 
action, the situation for the smoothing cycle is 
intermediate. Nevertheless, it should be pointed 
out that the global minimisation involved in in- 
verse blocking, at least naively, has its effect felt 
over the entire lattice. 

As evidence in favour of preserving long range 
physics, it was shown that the string tension 
is conserved under this smoothing (unlike for the 
method of ref. |Q, that reports changes up to 
25%). That smoothing is successful in removing 
noise is seen in fig. 4. Somewhat surprisingly close 
I-A pairs remained stable for distances as close as 
80% of the sum of their radii. Virtually no change 
is seen for up to 9 smoothing cycles. As I-A pairs 
are not solutions one would have expected some 
change. Possibly this is due to critical slowing 
down as can also occur for cooling (lij . It depends 
quite intricately on the details of the mapping 
implied by the smoothing cycle. 
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Figure 5. The SU{2) integrated instanton size 
distribution (integrated over bins of size 0.05 fm) 
for (3 = 1.4 (octagons), 1.5 (diamonds) and 1.6 
(squares) on a 8'* lattice, with resp. a — 0.188, 
0.144 and a==0.116 fm. From ref. [47]. 

In fig. 5 the resulting instanton size distribu- 
tion is given and we refer to ref. |47| for the 
details on the analysis. The instanton density 
is approximately 2.0 fm^^ and the size peaks at 
0.2 fm. This points to a rather dilute situation, 
although a study of correlations among the in- 
stantons seems to point to clustering. The sus- 



ceptibility xt = (230(10) MeV)4 agrees with the 
earlier value [Q. Cutting out instantons below 
0.27 fm (see fig. 3) would give Xt = (190 MeV)'', 
but a cut at half this value would cause no sig- 
nificant change. It would imply large scaling vi- 
olations for the improved cooling method, which 
were not observed in ref. [|o| . Both the suscepti- 
bility and the average size therefore disagree sig- 
nificantly with improved cooling. 

3.4. Spectral flow 

One can extract the topological charge by 
counting the number of chiral zero modes of the 
Dirac operator, using the Atiyah-Singer index 
theorem [|4|. Although only valid for smooth 
configurations, a lattice version could in principle 
be defined !§, Q = TOKpTr(75/(^TO))/7V/, with 
Kp a renormalisation constant that depends on 
the lattice definition of the Dirac operator Ip and 
75 used. Due to the discretisation of the lattice 
no exact zero modes exist. For Wilson fermions 
these would-be zero modes will typically be real 
and give rise to poles in the euclidean path inte- 
gral that are related to the "exceptional configu- 
rations" . Given the relation to instantons care is 
required in handling these configurations and in a 
recent series of papers a "pole shifting" algorithm 
has been proposed to remedy this problem ps[ |. 
One might also be tempted to extract the index 
by counting the real eigenvalues pO[ . 

Inspired by the overlap formulation for chi- 
ral fermions ]5l[ | a much simpler method to ex- 
tract the zero modes has been recently used |5^ . 
The Wilson Dirac operator has the property that 
H{U)=-f5{lp{U)-W{U) + 4-m) ishermitian. In 
the continuum this operator has a spectral flow 
as a function of m. Non-zero modes cross zero 
in pairs of opposite chirality (one going up, the 
other going down). A zero mode is generically 
responsible for an isolated crossing (the direction 
of crossing determines the chirality). It is these 
properties that are expected to be robust under 
discretisation . All that should change on the 
lattice is that the zero modes no longer cross at 
m = (we know that m needs to be tuned to 
rric ^ at flnite /3, but for individual configura- 
tions the crossings will occur at different values of 
to). The topological charge is now simply defined 
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as the number of net crossings. 
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Figure 6. Comparison of the SU{2) topolog- 
ical charge distribution from improved cooling 
(squares) and from the spectral flow (diamonds), 
on a 12"* lattice at /3 = 2.4. From ref. [52]. 

For very smooth instanton fields (instantons of 
large size), one should have near continuum be- 
haviour and the crossing of the zero mode should 
occur at small values of m. Perhaps the crossing 
value can thus be used to define the size of the 
instanton One way to study the correlation 
between size and crossing would be to use cooling 
to manipulate the size of an instanton. Essen- 
tial is that the spectral flow analysis is done on 
a coarse lattice. No smoothing is necessary. Fur- 
ther studies will be required, but the prospects 
are quite promising. 
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Figure 7. Comparison [52] for the fixed point ac- 
tion of the topological charge distribution using 
the spectral flow (Qievei) and after nine smooth- 
ing cycles (Qcycieg) on the same configurations. 

In fig. 6 a comparison is given for the topolog- 



ical charge distribution generated with the Wil- 
son action on the same lattice, measured with im- 
proved cooling [Q and with the spectral flow [^ . 
The agreement is excellent. From the spectral 
charge one extracts Xt = (184(6) MeV)^. 

Fig. 7 is based on 30 configurations generated 
with the fixed point action of ref. [^, on a 12"* 
lattice. The charge measured after nine smooth- 
ing cycles is compared with the spectral charge 
on the initial rough configuration [ |5^ . Cycling 
is seen to suppress the small charges. One con- 
figuration was used to trace the cause for this 
discrepancy [ ^ and it comes as an unpleasant 
surprise that the spectral charge was 1 without, 
3 after 9 and 2 after 12 smoothing cycles. 

4. EPILOGUE 

I report - you conclude. Who thought so much 
can be said about nothing. I humbly apologise 
to those that had hoped to find something else. 
I would have liked to discuss more on finite tem- 
perature and implications for the instanton liq- 
uid, on non-perturbative results in supersymmet- 
ric gauge theories as a testing ground for QCD 
and much more. But instantons are here to stay. 
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